Superradiant phase transitions with three-level systems 
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We determine the phase diagram of N identical three-level systems interacting with a single 
photonic mode in the thermodynamical limit (N — > oo) by accounting for the so-called diamagnetic 
term and the inequalities imposed by the Thomas-Reich-Kuhn (TRK) oscillator strength sum rule. 
The key role of transitions between excited levels and the occurrence of first-order phase transitions 
is discussed. We show that, in contrast to two-level systems, in the three-level case the TRK 
inequalities do not always prevent a superradiant phase transition in presence of a diamagnetic 
term. 



The collective superradiant coupling of a large number 
of atoms (or artificial atoms) has attracted a consider- 
able interest since the pioneering work of Dicke lj and is 
now the focus of many recent studies in cavity 2|-|6| and 
circuit quantum electrodynamics. The well-known 
Dicke model describes the coupling between a collection 
of two-level systems and a single photon mode. Remark- 
ably, for increasing light-matter coupling such a model 
predicts a superradiant phase transition with 
a doubly degenerate ground state above a critical vac- 
uum Rabi coupling. The so-called superradiant phase is 
characterized by a spontaneous polarization of the atoms 
and a spontaneous coherence of the cavity field in the 
ground state. In the case of time-independent Hamil- 
tonians, photons in the ground state are 'virtual' (i.e., 
bound in the cavity) and they cannot be radiated out of 
the cavity 1 1 3| unless some non-adiabatic modulation of 
the Hamiltonian is applied 14j (in analogy with the dy- 
namical Casimir effect [lol] ) . In this sense, the term 'su- 
perradiant' currently used for the Hepp-Lieb[l(| ground 
state is unfortunate, because it was originally introduced 
by Dicke [H for collective excited state radiative decay and 
not for ground state properties. Even in spite of no ex- 
tracavity emission from the ground state, occurrence of a 
superradiant critical point can be in principle monitored 
by measuring the dispersion of the collective excitations 
via standard optical techniques (e.g., through transmis- 
sion spectroscopy) . In the case of non-equilibrium phase 
transitions for pumped open systems, there is not a true 
ground state |2rl6|, but instead a stationary state, which 
is accompanied by the emission of real photons. 

However, for the case of electric dipole transitions, the 
Dicke model does not include the so-called diamagnetic 
term, which is proportional to the squared electromag- 
netic vector potential present in the minimal coupling 
Hamiltonian describing light-matter interaction in the 
non-relativistic regime. It is known that in the case of 
two-level real atoms, such a diamagnetic term is cru- 
cial, because it forbids the phase transition as a result 
of the TRK oscillator strength sum rule[8lll6j|. We point 
out that such no-go theorem cannot be necessarily ap- 
plied to time-dependent Hamiltonians with applied dress- 
ing fields [H, to the case of magnetic dipole coupling fl7j| 



and to more complex effective systems simulating the 
Dicke Hamiltonian with different degrees of freedom]^, 0] ■ 
Moreover, the no-go theorem^, [l6| is formulated for two- 
level systems, while the study of the multilevel case has 
been initiated only recently [l8|-|20| 

A generalized model including three-level atoms in 
the lambda configuration coupled to two photon modes 
has been theoretically investigated in a recent paper (2pj 
showing a very rich phase diagram with superradiant 
transitions of both second and first order (in the case of 
two-level atoms the superradiant transition is of second 
order). In such an interesting work[2(|, the diamagnetic 
term has not been included and the particular two-color 
lambda configuration has been considered, so the gen- 
eral case of an arbitrary three-level system needs to be 
explored. In a recent letterflij], it has been claimed that 
for multilevel atoms coupled to a single photon mode 
the no-go theorem can be generalized using again the 
TRK oscillator strength sum rule: in the proof reported 
in such a work 18} , through a perturbative argument, it 
is assumed that transitions between excited states can 
be always neglected in the thermodynamical limit: this 
is rather surprising [l9j, since in the work by Hayn et 
al.[20|, transitions between excited states play instead 
a crucial role in the thermodynamical limit and are re- 
sponsible for the first-order transition boundaries. It is 
therefore a fundamental problem to explore superradiant 
phase transitions with arbitrary three-level systems. 

In this letter, we investigate the existence of super- 
radiant phase transitions of a system consisting of N 
three-level systems coupled to a single photon mode in- 
cluding the diamagnetic term in the Hamiltonian and 
TRK inequalities. We show the rich phase diagram in 
the thermodynamical limit (N —¥ oo), obtained via a 
multilevel Holstein-Primakoff approach. We have stud- 
ied particular configurations (ladder, lambda, V-type) 
and also the general three-level coupling configuration. 
We find that for three-level systems TRK inequalities do 
not always prevent superradiant phase transitions and 
that excited state transitions play a crucial role in the 
thermodynamical limit. As sketched in Fig. 1, 

let us consider N identical three-level systems, whose 
states are {|0fc),|lfc),|2fc}} , (k = 1,2,..,N) where k is 
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FIG. 1: A sketch of the considered system consisting of N 
identical three-level atoms identically coupled to a single cav- 
ity mode. Each three-level system is represented by the tran- 
sition frequencies wio = wi — ujq > 0, CJ21 = U2 - > 0, 
oj2o = W21 + wio > and by the oscillator strengths /oi > 0, 
/12 > 0, fo2 > (see definition @ in the text). 



FIG. 3: V-type configuration (/12 = 0). Photonic order pa- 
rameter as a function of the oscillator strengths /01 and /02 • 
The red-dashed lines indicate the boundaries imposed by the 
TRK sum rule, namely < foi + /02 < 1. Parameters: 
D = ujcav, wio = iOcav, W21 = O.liOcav ■ In this configura- 
tion, the superradiant phase is always incompatible with the 
TRK boundaries. 
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FIG. 2: Results for the ladder configuration (/02 = 0). Pho- 
tonic order parameter x = a/\fN — (a)/y/N in the ground 
state as a function of the oscillator strengths /01 and /12 for 
each three-level atom. The normal phase (N) is character- 
ized by x = (black color). In the superradiant phase (SR), 
there is a spontaneous coherence i/O. Note that a discon- 
tinuous jump of x denotes a first-order superradiant phase 
transition, otherwise the transition is of second order. The 
area below the red-dashed lines indicate the region described 
by the TRK inequalities < /01 < 1 and < /10 + /12 < 1 
where /10 = —foi. Parameters: D — 3uj cav , u)io = 0.1u C ai>, 
CJ21 = ijJcav The collective vacuum Rabi frequency fioi an d 
Hi2 are obtained through the relationship in Eq. JSJ. 



the atomic index. Each atom is assumed to be indepen- 
dent from the others and to interact identically with a 
single bosonic, photonic mode. The energies of the three 
levels are fiwo < hui < hu>2- 

By introducing the collective operators Sy = 



FIG. 4: Results for the lambda configuration (/01 = 0). Pho- 
tonic order parameter versus /12 and /o2- The area below 
the red-dashed lines is compatible with the TRK inequalities, 
namely < /12 < 1 and < /02 < 1- Parameters: D = 3u} ca v, 
cjio = Q.lui C av, W21 = 0-9uj cav . As in the ladder case, the su- 
perradiant part of the diagram has an overlap with the region 
compatible with the TRK inequalities. 



SfcLi \ik)(jk\, the light-matter Hamiltonian reads 

D{a 



at) 2 



(1) 



The cavity mode is described by the frequency u) cav and 
by its creation (annihilation) bosonic operator a' (a). 
The coupling between the cavity mode and the atomic 
i — > j transition is quantified by the collective vacuum 
Rabi frequency Sly. Qij/^/N is the light- matter coupling 
per atom. The term proportional to D and quadratic in 
the photon operators is the so-called diamagnetic term. 
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FIG. 5: Top panel: the filled red volume in the (/oi, fo2, /12) 
space represents the superradiant part of the diagram com- 
patible with the TRK inequalities < foi + /02 < 1 and 
< fio + fi2 < 1. Parameters: D = 5ui cav , wio = 0.17(jj cav , 
U21 — Wca»- Middle panel: the photonic order parameter x 
as a function of /01 and /02 for a fixed value /12 = 0.735. 
The area below the red-dashed line corresponds to the re- 
gion satisfying the TRK inequalities. Bottom panel: the 
black line represents the spatial profile of an illustrative arti- 
ficial one-dimensional potential with squared wells (as those 
realizable with semiconductor heterostructures) in units of 
the energy E c 



1,- 



, ft where the spatial coordinate is ex- 
pressed in units of the length L. The horizontal lines de- 
pict the energies of the three bound states with their corre- 
sponding wavefunctions (other states are in the continuum). 
For this potential shape, one obtains the oscillator strengths 
(,foi = 0.3995, f 02 = 0.40 69, f 12 = 0.735) and anharmonicity 
ratio W10/W21 = 0.1709 corresponding to the point depicted 
by the blue cross in the middle panel. 



electron-light interaction in the so-called p • A gauge, 
where p is the electron momentum operator, q the elec- 
tron charge and A is the electromagnetic vector potential 
operator. As it can be deduced by the treatment in Ref. 
[8|, the value of the vacuum Rabi frequencies fiy are 
linked to the diamagnetic term amplitude via the oscil- 
lator strengths as follows: 



— = fijD 



j< 



where w,- 



— u)i is the transition frequency and 
fij = _ Ui)\dij\ 2 



(2) 



(3) 



is the transition oscillator strength with the corre- 
sponding atomic electric dipole matrix element (along 
the direction of the photon mode polarization). Indeed, 
as shown in Ref. [1] , HD = ^n e iNAl where q is the 
electron charge, n e i the number of electrons per atom and 
A = Ao(a + at) is the electromagnetic vector potential in 
the position where the atoms are assumed to be coupled 
(the spatial variation of the cavity field is neglected in 
the region occupied by the atoms, i.e., the electric dipole 
approximation has been considered). For a given col- 
lection of atoms and for a given cavity, D is a constant 
depending on the density of atoms in the cavity volume 
via the factor AAq. A related similar model can be ob- 
tained by considering a two-dimensional electron gas in a 
semiconductor quantum well heterostructure 2l|-|24| (re- 
placing the bare electron mass with the semiconductor 
conduction band effective mass). 

The TRK oscillator strength reads J2j fij = 1- Note 
that by definition the oscillator strengths /y > if 
ujj > u>i, while it is negative if u)j < Wj. Hence, for a tran- 
sition from the ground level |0) to the first excited level 
|1), we have always < /01 < 1 as foj > 0. For two-level 
system models, the presence of the diamagnetic term and 
the constraint < /01 < 1 is sufficient to prevent the su- 
perradiant phase transition (no-go theorem) [l6| . Here, 
we will analyze what happens with three-level systems. 
By applying the method detailed in Ref. [25[ , we can 
express the generalized collective transition operators by 
using a multilevel Holstein-Primakoff transformation in 
terms of bosonic annihilation (creation) operators Ck ( 
cj,). When we use this procedure for the three- level case, 
we have to choose a state of reference that we will call r. 
The multilevel boson mapping reads: 



Sfefe = c\c k (k ^ r) 



V k ^r 



E rr = n - ^2 4a 



(4) 
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The collective transition operators £y defined above do Note that /?o+/?2 — N (which is equivalent to y 2 +z 2 < 1) 



not commute and are such that 



2-Hj , S 



5iiY.kj. The Hamiltonian can be rewritten as 
H/H =uj cav a)a + ^^(wfe — (jj r )c\.Ck + co r N 

k^r 



ik^r 



Cl c k 



i>j 



(a + a f ) + D(a + a f ) 2 



(5) 



In the following, we choose r — 1 as reference state 
(choosing another reference state leads to the same re- 
sults). In order to determine the phase diagram in the 
thermodynamical limit, we can use the mean-field ap- 
proach as in Refs. [lj, by the replacement (a) = 
(a 1 *) = a , (cj) = (cj) = fij with j 6 {0,2}. It can be 
shown that due to the form of Eq. ([!} the solutions for 
a and /3j are necessarily real numbers. Therefore, we ob- 
tain the mean-field expression of the ground state energy 

© 

E G /h ={u cav + AD)a 2 + lo 01 PI + UJ2102 + ^iN 



Aa 



(6) 



whose global minimization will give the values of a (pho- 
tonic coherence), /?o and ft (collective atomic coher- 
ences) in the ground state. A photonic order parameter 
a = implies that the system in the Normal (N) phase. If 
a / 0, a SuperRadiant (SR) phase occurs [12|. In order 
to minimize the ground state energy in the thermody- 
namical limit [N — > oo) it is convenient to introduce the 
rescaled quantities x = , y — and z = The 
ground state energy can therefore be rewritten as 

E G /h 



N 



--{u cav + AD)x +u} 01 y +uj 2 iz +u>! 



Ax 



(Qioy + n 2 iz)vi -y 2 - z 2 + tt 2 oyz 



As 



2 E G 



2(w c 



(7) 

AD) > 0, the minimization with 
respect to x is therefore trivial, giving a; as a function of 
y and z. In order to minimize the ground state energy, 
we therefore need to minimize with respect to y and z 
the following energy function (obtained from Eq. (|7P): 

£g =uoiy 2 + U21Z 2 + wi 
4 



AD 



(fiioj/ + fi 2 i 2) V 1 - y 2 - + ^202/2 

(8) 



must be fulfilled within the Holstein-Primakoff frame- 
work. Hence, the global minimum of £q has to be looked 
for in a disk with unity radius centered in the origin of 
the Oyz plane, which can be done with a straightforward 
numerical calculation. 

IWe start by considering the ladder configuration (/02 = 
0) for the three-level systems. In Fig. [51 we have plotted 
results for the photonic coherence x = a/VN = (a)/y/N 
as a function of the oscillator strengths /01 and /12 of 
the ladder coupling for D = 3uj cav , u>2i — uj cav and 
W10/W21 =0.1. We point out that these parameters cor- 
respond to a strong anharmonicity with the photon mode 
in resonance with the excited transition. The value of the 
collective vacuum Rabi frequencies is given by Eq. @. 
The normal phase occurs when a = (black region). 
The results in Fig. [2] show that in the considered sys- 
tem a superradiant phase (a ^ 0) do occur. The phase 
transition boundary can be of both of first and second 
order. The first-order phase transition occurs when a 
has a discontinuous jump from to a finite value (upper- 
part of the frontier) . The first-order transition boundary 
is due to the excited transition 1 — > 2 as in the case re- 
cently studied by Hayn et al. 20] : at the transition there 
is a macroscopic occupation of the intermediate state 
The second-order phase transition occurs when a is not 
discontinuous, but its gradient is. The important point 
to consider here is the fact that due to the TRK sum 
rule the oscillator strengths are subject to constraints. 
For the ladder configuration, we have < /01 < 1 and 
< /10 + /i2 < 1, where f w = -f i < 0. In Fig. [5] the 
TRK boundaries imposed by such inequalities have been 
indicated: the area compatible with the TRK sum rule is 
below the red-dashed line. It is therefore apparent that 
there is an overlap between the superradiant part of the 
diagram and the region compatible with the TRK sum 
rule. Note that the relevant superradiant phase region 
compatible with the TRK inequalities contains a first- 
order transition boundary. 

In Fig. [21 we consider instead the V-type configuration, 
where /12 = 0, i.e., there is no coupling between the ex- 
cited states. In such a configuration, the photonic order 
parameter x = a/^/N is plotted as a function of /01 and 
/02 (parameters in the caption). The TRK sum rule im- 
poses the inequality < /oi + /o 2 < 1, again indicated by 
the red-dashed line. In such a V-type configuration, there 
is a superradiant phase boundary, which is always of sec- 
ond order. In fact, the transition between excited states 
is by definition inactive in such a configuration. Impor- 
tantly, we notice that here there is no overlap between 
the superradiant part of the diagram and the area com- 
patible with the TRK sum rule. Indeed, the TRK sum 
rule always prevents the superradiant phase transition in 
the V-type configuration as in the two- level caseQ. 

In Fig. m we show results for the lambda configura- 
tion (/01 = 0): the photonic coherence is shown versus 
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/i2 and fo2 (parameters in the caption). Here, the TRK 
inequalities < /12 < 1 and < /02 < 1, whose bound- 
aries are once again delimited by the red-dashed line. As 
in the ladder case, in the lambda configuration there is 
a superradiant phase in the region compatible with the 
TRK inequalities (are below the red-dashed line) with 
the transition boundaries being of the first order. 

So far, we have shown the simplest three-level config- 
urations (ladder, V-type and lambda). In Fig. El we 
show results for a generic three level system where all 
the three oscillator strengths are finite (detailed param- 
eters in the caption). The superradiant part satisfying 
the TRK inequalities is depicted by the red filled volume 
in the (/01, /02, /12) space (top panel). The middle panel 
shows the photonic order parameter on a planar section 
of such oscillator strength three-dimensional space. We 
point out that in the generic three-level case the interest- 
ing superradiant part compatible with the TRK inequal- 
ities is widened due to the additional freedom associated 
to the third oscillator strengths . In the bottom panel of 
Fig. El we give an illustrative example of a spatial arti- 
ficial potential providing a three-level system with oscil- 
lator strengths (/01, /02, /12) and anharmonic spectrum 
corresponding to the cross in the middle panel, for which 
a superradiant transition is possible while satisfying the 
TRK constraints. As a perspective, it will be interesting 
in the future to address more complex multilevel struc- 
tures and to investigate also the role of direct Coulomb 
interactions [26j . 

In conclusion, we have determined the phase diagram 
for a model system consisting of N three-level systems 
coupled to a single photonic boson mode, by including 
the diamagnetic contribution in the light-matter cou- 
pling. We have demonstrated that in the considered 
model system superradiant phase transitions can occur 
while preserving the TRK inequalities, in stark contrast 
to the case of two-level systems. We have found that 
the transition between excited levels have a key role in 
such superradiant phase transition in contrast to what 
assumed in Ref. [F8[. Our results show that the physics 
of superradiant phase transitions with multilevel systems 
might be achievable in a broad range of physical systems, 
especially those where it is possible to engineer the spec- 
tra and oscillator strengths. We would like to thank T. 
Brandes, I. Carusotto, C. Emary and J. Keeling for stim- 
ulating discussions. C. C. is member of Institut Univer- 
sitaire de France (IUF). We acknowledge support from 



the ANR project QPOL. 
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